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We have recently developed a quantized fluctuational electrodynamics (QFED) formalism to de¬ 
scribe the quantum aspects of local thermal balance formation and to formulate the electromagnetic 
field ladder operators so that they no longer exhibit the anomalies reported for resonant structures. 

Here we show how the QFED can be used to resolve between the left and right propagating fields 
to bridge the QFED and the quantum optical input-output relations commonly used to describe 
selected quantum aspects of resonators. The generalized model introduces a density of states con¬ 
cept describing interference effects, which is instrumental in allowing an unambiguous separation 
of the fields and related quantum operators into left and right propagating parts. In addition to 
providing insight on the quantum treatment of interference, our results also provide the conclusive 
resolution of the long-standing enigma of the anomalous commutation relations of partially confined 
propagating fields. 


I. INTRODUCTION 

It has recently been suggested m that in contrast 
to earlier predictions the commutation relations of pho¬ 
ton ladder operators have directly measurable physical 
significance. This is especially interesting in the case 
of resonant structures where conventional theoretical de¬ 
scriptions have been shown to lead to anomalous com¬ 
mutation relations of the general form [a,d^] = A yf: 1, 
instead of the canonical form with A = 1 The 

anomaly is a direct consequence of the conventional piece- 
wise normalization of the optical modes in respective ho¬ 
mogeneous regions and of the interference effects coupling 
the modes propagating in different directions. Originally 
these anomalies were argued to bear no physical signifi¬ 
cance, but recently it has been shown that the anomalous 
commutation relations lead to the existence of a thresh¬ 
old for second harmonic generation when it occurs inside 
microcavities mm- In addition, anomalies in the com¬ 
mutation relations have also been shown to prevent sys¬ 
tematic description of the local thermal balance between 
the field and interacting media Eiiiiin]. Experimental 
measurements of the onset of the second harmonic gen¬ 
eration or the thermal balance formation in microcav¬ 
ities could therefore confirm the theoretical predictions 
that the conventional mode normalization introducing 
these anomalies is not sufficient for field quantization of 
resonant structures. Here, we develop a field quantiza¬ 
tion approach that (1) fully eliminates the anomalies for 
propagating fields, (2) bridges the classical propagating 
wave and commutation-relation-preserving quantum de¬ 
scriptions, and (3) allows formulating conceptually sim¬ 
ple models for optical energy transfer and the formation 
of thermal balance in interfering nanostructures. 

One of the most widely used quantization approaches 
for describing spatial field evolution in resonant struc¬ 


tures is the input-output relation formalism of the pho¬ 
ton creation and annihilation operators. The formalism 
was originally developed for dispersionless and lossless 
media m and later extended for lossy and dispersive 
dielectrics by several groups [Mni- The quantization 
procedures studied, e.g., by Barnett et al. [H] clearly 
highlight that the noise and field operators in nonuni¬ 
form systems are position dependent and that the vector 
potential and electric-field operators obey the well-known 
canonical commutation relation as expected [HUH]. The 
canonical commutation relations in these early models 
did not, however, extend to the ladder operators which 
were found to exhibit anomalies in resonant structures 
[5|. The anomalous commutation relations of the ladder 
operators were later studied in several works [SHni but 
no clear resolution for the anomalies was found. Instead 
it was concluded that the anomalies as well as the exact 
form of the ladder operators within resonant structures 
were irrelevant as long as the field commutation relations 
and classical field quantities were well defined. However, 
this made it impossible to fully quantize the optical fields 
in resonant structures. 

To shed more light on the anomalous commutation 
relations, we have very recently developed a quantized 
fluctuational electrodynamics (QFED) scheme based on 
generalizing the fluctuational electrodynamics to quan¬ 
tum optical fields laiiiHi. Using the QFED approach 
we were able to formulate the canonical commutation 
relations preserving ladder and photon-number opera¬ 
tors for the total electromagnetic (EM) field [3 H] [TU] . 
However, even in the QFED framework, it has not been 
evident how to separate the ladder and photon-number 
operators to left and right propagating parts, which is 
also essential for the final resolution of the anomalies 
and for bridging the classical propagating wave descrip¬ 
tions and the commutation-relation-preserving quantum 
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descriptions. In this work, we show that the QFED can 
be extended to resolve between the left and right prop¬ 
agating fields, fully preserving the canonical commuta¬ 
tions with A = 1 also for the left and right propagating 
field ladder operators. As it turns out, the separation to 
the left and right propagating helds becomes possible and 
conceptually simple when one introduces a new density 
of states concept describing the fundamentally important 
interference effects. The added insight obtained using 
the introduced concepts provides a more detailed under¬ 
standing of the quantization of optical fields in complex 
surroundings, and can be used, e.g., for finding simple 
photon-number-based expressions for the quantum opti¬ 
cal Poynting vector as well as more detailed description 
of quantized energy flow in resonant structures. 


II. FIELD QUANTIZATION 
A. Photon operators 


In contrast to previous approaches, the requirement of 
the QFED is the preservation of the local canonical com¬ 
mutation relation [a{x,uj),a^{x^uj')] = i5(a; — uj') also in 
resonant media nun]. This requirement leads to con¬ 
ceptually simple definitions for the position-dependent 
ladder and photon-number operators as a weighted sum 
over the incident fields and the noise. For the expecta¬ 
tion value of the photon-number operator, for instance, 
this weighted sum reads as [2] 


{n{x,U})) 


jToa PNL{x,uj,x'){fi{x',uj))dx' 
j!^^PNL{x,UJ,x')dx' 


( 1 ) 


where {ri(x,ui)) is the source field photon-number expec¬ 
tation value which for thermal helds is given by the Bose- 
Einstein distribution {ri{x',u})) = — 1] 

with the position-dependent temperature of the medium 
given by T{x'). The weighting coefficient pNL(a;, w, cc') in 
Eq. 0 is given by 


PNL(a;,w,x') 


^'^\£{.x,w)\ 

TTd^S 


£i{x',Lu) 


X {|G(a:,w,x')p -I- 


dG{x,io,x') 
k{x, uj)dx 



where c is the speed of light in vacuum, S is the area of 
quantization in the y-z plane, e{x,uj) = n(x,a;)^ is the 
relative electric permittivity of a nonmagnetic medium 
with refractive index n(x,uj) and ei(x,a;) is its imagi¬ 
nary part, k(x,ui) = uin{x,u!)fc is the wave number, and 
G{x, UJ, x') is the Green’s function of the Helmholtz equa¬ 
tion given for selected layered structures in Ref. |3l The 
quantity p^i,(x,uj,x') is here referred to as the nonlocal 
density of states (NLDOS) since it highlights the nonlo¬ 
cal origin of the local density of states (LDOS) p{x,uj) 


given as 


p{x,uj)= pj<ii^{x,uj,x')dx', (3) 

J —OO 

and appearing in the denominator of Eq. 0. After inte¬ 
gration, the LDOS can also be expressed in the more 
familiar form in terms of the imaginary part of the 
Green’s function [4]. Using this definition, the NLDOS 
PNL(a^)W,x') accounts for both the electric-field (term 
|Gp) and the magnetic-field (term \dG/{kdx)\^) contri¬ 
butions. 


B. Quantum optical Poynting vector 


The quantum optical Poynting vector S{x,t) that will 
be used as a starting point for separating the field 
components is defined in terms of the positive (-I-) 
and negative (—) frequency parts of the electric- and 
magnetic-field operators E{x,t) and B{x,t) as S{x,t) = 
eoc^[E~{x,t)B~^{x,t) + B~{x,t)E~^{x,t)] [HI[19]. Using 
the QFED framework and substituting the electric- and 
magnetic-field operators as given in Ref. [3| we write the 
Poynting vector expectation value at angular frequency 
UJ as 

/ OO 

Pif(x, uj, x'){f}{x',uj))dx', 

-OO 

(4) 


where v(x,uj) = c/nr(x,w) is the energy propagation ve¬ 
locity, nr{x,uj) is the real part of the refractive index 
n(x,uj) = y^e(x,uj), and 


, 2ujnr{x,uj) , 
Pif(x,uj,x ) = ——£i(x ,UJ) 


X Re(^iujG{x,uj,x') 


, 9G* (x,uj,x') 


dx 


)• (5) 


The quantity pif(x,uj,x') at held point x essentially de¬ 
scribes the contributions of the left and right propa¬ 
gating helds originating from the source point x'. The 
term pif{x, uj, x') fully accounts for the rehections, losses, 
and interference and is closely related to the concepts of 
LDOS and NLDOS. Therefore, we refer to it as the in¬ 
terference density of states (IFDOS). In contrast to the 
NLDOS, the integral of the IFDOS with respect to x' 
is always zero as required, e.g., by the fact that in a 
medium in equilibrium, there is no net power how, i.e., 
{S{x,t))i^ = 0 when {fi{x',uj)) is constant. 


C. Left and right propagating fields 

To generalize the LDOS and photon-number concepts 
in Refs. |3| and |2| and Eq. 0 to separately account 
for the left and right propagating fields, we write the 
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left and right propagating field Poynting vector expecta¬ 
tion values {S+{x,uj)) and {S-{x,uj)) as {S±{x,uj))ui = 
tiujv{x,uj)p±{x,uj)(^{n±{x,uj)) + i), where p±{x,uj) and 
{n±{x, oj)) are the left and right propagating field LDOSs 
and photon numbers to be determined, and the term one 
half describes the zero-point fluctuation current. The left 
and right propagating photon numbers must additionally 
satisfy two equations: the total Poynting vector must be 
given by 


{S{x,t))^ = Hujv{x,oj)p+{x,uj)(^{n+{x,oj)) + 

- fujv{x,uj)p-{x,uj)(^{n_{x,uj)) + 

( 6 ) 

and the total energy density {u{x,t))^ = 

hujp{x,u})({n{x,uj)) -b 1) m by 

{u{x,t))^ = hujp+{x,uj)(^{n+{x,uj)) + 

+ fiwp-{x,u})(^{n_{x,oj)) + ^y ( 7 ) 

At zero temperature, where {n^{x,ui)) = {h-{x,uj)) = 0, 
the Poynting vector is zero and thus p+{x,uj) = p-{x,ui) 
in Eq. (§. Respectively, Eq. 0 at zero temperature 
leads to the relation p^(x,uj) + p-{x,ui) = p{x,uj). To¬ 
gether, these conditions uniquely define the left and 
right propagating LDOSs in terms of the total LDOS 
as p+{x,uj) = p-{x,uj) = p{x,L>j)l2. 

Using the above local density of states relations, we 
can uniquely solve the left and right propagating pho¬ 
ton numbers from Eqs. Q and Q as {n±{x,u})) = 
[hujp{x,uj)]~''-[{u{x,t))^ ± {S{x,t))^/v{x,u})] - 1/2. In 
terms of the source field photon number this corresponds 
to 


{n±{x,uj)) 


f^^[PNh{x,uj,x') ± piF{x,uj,x')]{fi{x\uj))dx' 
J^^[PNLix,UJ,x') ± piF{x,UJ,x')]dx' 


Equation ^ shows that the propagating field photon- 
number expectation values are also obtained as a 
weighted sum of the source field values, but the weight 
factor now includes an additional term describing the in¬ 
terference and propagation direction. In the denomina¬ 
tor, one can neglect pi-p{x,ui,x') as it integrates to zero 
indicating that the denominator is simply equal to the 


LDOS in Eq. 


Above we have only focused on the photon-number ex¬ 
pectation values that can be directly extracted from the 
Poynting vector. To find the corresponding ladder and 
photon-number operators in the QFED we will further 
investigate the forms of the photon annihilation opera¬ 
tors d+(x,uj) and d-(x,uj) that lead to the expectation 
values in Eq. 0 and fulfill the canonical commutation 
relations. The photon annihilation operators fulfilling 


these conditions are of the form 

1 i 


d±{x,u) = 


i(^±7r/4) 


y/pix,u;) 7-00 

X ^/pf^F{x,UJ,x') ± pif(x,w,x') f{x',uj)dx', 

(9) 


where f{x',ui) is a bosonic source field opera¬ 
tor obeying the canonical commutation relation 
[f{x,oj),p{x',oj')] = S{x — x')5{uj — uj') and which 
is related to the source field photon number as 
{f){x',uj)) = J{p{x',uj)f{x",uj'))dx"duj' [3J|31[Tn]. The 
phase factor is in principle arbitrary and it does not 
play a role in our calculations as it cancels in the com¬ 
mutators. The total field annihilation operator d{x, u) 
is given by the sum d{x,uj) = -^[d+{x,uj) + d_(a;,a;)]. 
It is straightforward to check that the left and right 
propagating field annihilation operators in Eq. 0 
also obey the commutation relation of the form 
[a±(a;,a;),a^(a:,w')] = 5{uj — uj'). With these choices, 
however, the cross-commutators become nonzero as 
[a±(a;,w),d^(a;,w')] ^ 0 due to the coupling of the 
left and right propagating fields originating from the 
same source points by the reflecting interfaces. This 
cross-commutator form is intuitively reasonable and 
does not appear to present any complications as the only 
commutation relations directly linked to the studied 
physical observables are the self-commutators. 

As shown by Eqs. ([^ and 0 it is necessary to sep¬ 
arately account for all the individual source points and 
their mutual interference to arrive to the correctly com¬ 
mutating operator forms. Similar book keeping is also 
present in classical fluctuational electrodynamics (FED). 
In contrast to the FED, however, to describe the quan¬ 
tum features the photon ladder and number operators 
need to be renormalized to fully satisfy the commutation 
relations. 

In addition to describing the total energy density and 
energy flow presented in Eqs. 0,0, and 0, the QFED 
formalism is also capable of separating the total field pho¬ 
ton numbers to their local electric- and magnetic-field 
equivalents {ne(x,uj)) and {h^{x,uj)) that are responsi¬ 
ble for direct interactions with materials and determine, 
e.g., the self-consistent local temperature of the interact¬ 
ing media as discussed in Refs. 0] and [31 Essentially these 
electric- or magnetic-field specific quantities and the cor¬ 
responding LDOSs can be obtained by using Eqs. 0 
and 0 when only the electric- or magnetic-field term 
in the NLDOS in Eq. 0 is taken into account 0]. The 
electric- and magnetic-field specific quantities were previ¬ 
ously shown to have quite distinct properties as compared 
to the total field quantities and to include, e.g., oscilla¬ 
tions in the field temperatures 011 ]. In the propagating 
operator formalism, however, the electric- and magnetic- 
field specific ladder operators are again united with the 
total propagating field operators because the direct in¬ 
terference effects between the left and right propagating 
fields have been eliminated when projecting the ladder 
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FIG. 1. (Color online) Optical cavity consisting of three ho¬ 
mogeneous layers. We calculate the left and right propagating 
field photon-number expectation values in each layer. 



operators to the left and right propagating operators. 
This fully agrees with our previous results HE]: in the 
present formalism the separation to left and right prop¬ 
agating fields also fully separates the interference effects 
from the local fields, whereas the formalism simultane¬ 
ously capturing both left and right propagating fields in 
a single term must also capture the interference effects. 


FIG. 2. (Color online) Left and right propagating photon 
numbers (fi-) and (h+) and the total photon number (h) in 
a lossy cavity structure at the first resonant energy ftta = 
0.046 eV (A = 26.9 ^m). The media from left to right have 
refractive indices ^eT = 2.5 -I- OAi, = 1.2 -|- 0.2i, and 
= 1.5 + 0.5i, and source field temperatures Ti = 300 K, 
T 2 = 200 K, and Tg = 100 K. 


III. RESULTS 

To better illustrate the physical implications of the pre¬ 
sented concepts we briefly discuss the properties of pho¬ 
ton numbers of the left and right propagating fields and 
compare them to the corresponding total field photon 
number in an optical cavity consisting of three homoge¬ 
neous layers as illustrated in Fig. 


can substantially increase the field magnitude inside a 
resonator, the photon-number values inside the cavity 
and at the outputs in Eq. (10 1 are always between the 
input field photon numbers. This also ensures that in 
global thermal equilibrium all the photon numbers are 
equal and no photon-number accumulation can occur in¬ 
side the cavity. 


B. Lossy cavity structure 


A. Lossless cavity structure 

In a lossless configuration, the left and right propagat¬ 
ing field photon numbers are piecewise continuous and 
only depend on the cavity geometry and the input fields 
(fii+) and (u-a-) incident from the left and right. In dif¬ 
ferent regions, they can be written as 

(ni_) = |7^lp(nl+) -b x/iyTea |Ti'7^'P(n3-), 

^ \/£ 2 /gi |Tip(hi+) -b ^ 82/63 |r27^;p(n3-} 

“ Re[l -b 

\ \/£2/£i |Tl7^2p(nl+) + a/£2/£3 IT 2 P(n3-) 

~ Re[l-b 27^i7^2^^2e2*'==‘^2] 

(n3+) = \Je-ilei |T(7^P(hi+) -b |7^2p(n3-), ( 10 ) 

where ^2 is the cavity thickness, is the wave num¬ 
ber inside the cavity, V 2 = 1/(1 + 7?.i = 

(ri -b r 2 e^'‘^^'^^)v 2 , 7^2 = ra, Ti = tiV 2 , T 2 = ^ 2 , = f'l, 

^2 = (^2 + '^ 2 )^ 2 , T( = t'l, and T-i = ^ 2^2 with the 

conventional single interface Fresnel reflection and trans¬ 
mission coefficients for left incidence and U, i G { 1 , 2 }, 
and right incidence r' and t', i S {1,2}. In contrast, 
e.g., to the electric-field values where resonance effects 


In a lossy structure the photon numbers are no longer 
piecewise constant and all material points can act as field 
sources through the source field (r)(x',w)), which is re¬ 
lated to material temperature. To illustrate this, we 
study a lossy cavity structure, where the refractive in¬ 
dices of the media are y/^ = 2.5-b0.4z, y/^ = 1.2-b0.2i, 
and y/^ = 1.5 -b 0.5*, and the layer temperatures are 
Ti = 300 K, T 2 = 200 K, and T 3 = 100 K. Here the 
layer temperatures are set to constant predefined values 
for simplicity even if the QFED formalism also allows 
calculating the in-cavity temperature distribution self- 
consistently if the studied layers were to be considered as 
thermal insulators Ellin]. Figureshows the the total, 
right propagating, and left propagating photon numbers 
(n), (h_|_), and (n_) as a function of position at the first 
cavity resonance tiiu = 0.046 eV (A = 26.9 i*m), where 
the layer temperatures correspond to steady-state pho¬ 
ton numbers 0.20, 0.074, and 0.0048. The photon num¬ 
bers are highest at the leftmost medium at Ti = 300 K 
and decrease towards the rightmost medium at T'i = 100 
K. The right propagating photon number notably de¬ 
creases at and after the first interface due to reflection 
and thermalization, eventually reaching equilibrium with 
the lossy medium in the rightmost layer. The left propa- 
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FIG. 3. (Color online) (a) Total LDOS, (b) Tes, (c) and (d) in a lossy cavity structure. The media from left to right 
have refractive indices y/ei = 2.5 + 0.4i, = 1.2 + 0 . 2 z, and ^/e 3 = 1.5 + 0.5i, and layer temperatures Ti = 300 K, T 2 = 200 

K, and T 3 = 100 K. Solid lines denote the boundaries of the cavity and dashed lines denote resonant energies. The LDOS is 
given in the units of 2/{'kcS). 


gating photon number notably changes at the interfaces 
and in the middle and leftmost layers. It can be also 
clearly seen that the total photon number is the average 
of the left and right propagating photon numbers as ex¬ 
pected, since the photon number essentially describes the 
average photon number in the collection of optical modes 
under study. 

As the photon-number expectation value depends 
strongly on the frequency, it is convenient to illustrate 
the results by using the effective field temperature that is 
defined in terms of the photon-number expectation value 
a,s Ted{x,uj) = Hu}/(^kB\n[l + l/{h{x,uj))]) [U[Tn]. This 
corresponds to the steady-state temperature of a small 
temperature probe interacting only with a single mode 
d]. Figure shows the total LDOS and field tempera¬ 
tures corresponding to the total, right propagating, and 
left propagating fields as a function of position and pho¬ 
ton energy. In contrast to the case of a lossless structure, 
the field quantities are position dependent. The total 
LDOS in Fig. [^a) is also oscillatory inside the cavity 
and reaches its maxima at resonant energies hw = 0.046 
eV (A = 26.9 /rm), huj = 0.097 eV (A = 12.7 i^m), and 
hoj = 0.150 eV (A = 8.29 /xm). Also in the left- and 
rightmost layers, the total LDOS is position dependent 


and oscillatory near interfaces. The oscillations of the 
LDOS follow from the interference effects combined with 
the material polarizability in analogy with the Purcell 
effect [4]. 

Despite the oscillations in the LDOS, the total effective 
field temperature Teg in Fig. |^b) and the effective field 
temperatures TJj and corresponding to the right and 
left propagating fields in Figs.j^c) and[^d) decrease to¬ 
wards the right medium at lower temperature similar to 
the photon numbers in Fig. On the left and right Tgg, 
TgJj, and also asymptotically approach equilibrium 
values corresponding to material temperatures. When 
compared to T^, the magnitude of is everywhere 
lower since the source field temperature on the right is 
lower than the source field temperature on the left. 


IV. CONCLUSIONS 

In conclusion, we have developed a generalized quan¬ 
tum optical noise formalism QFED that can unambigu¬ 
ously describe the quantum aspects of propagating opti¬ 
cal fields in arbitrary stratified media, while being fully 
compliant with the canonical commutation relations. 
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In particular, the QFED allows calculating position- 
dependent photon-number expectation values for the left 
and right propagating helds and fully eliminates the 
anomalies of the ladder operators in optical cavities. 
In our model, the commutation relations are therefore 
canonical. This implies that, in contrast to previous mod¬ 
els that involve anomalous commutation relations, our re¬ 
sults do not predict, e.g., any observable threshold for the 
second harmonic generation inside cavities as the thresh¬ 
old is directly linked to the ladder operators. Experimen¬ 
tal measurements of the existence of a second harmonic 
generation threshold may therefore allow demonstrating 
the importance of correct normalization of the commuta¬ 
tion relations as well as the affiliated normal modes. In 
addition, the QFED framework enables, e.g., the separa¬ 
tion of the quantum optical Poynting vector and related 


field quantities to their left and right propagating com¬ 
ponents using a photon-number-based presentation and 
the interference density of states. In practical modeling 
tasks, the QFED provides simple tools for studying opti¬ 
cal energy transfer and the formation of thermal balance 
in complex interfering nanostructures, highlighting the 
fundamentally nonlocal nature of the energy transfer. 
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